Finite temperature stability and dimensional crossover of exotic superfluidity in 

lattices 
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We investigate exotic paired states of spin-imbalanced Fermi gases in anisotropic lattices, tuning 
the dimension between one and three. We calculate the finite temperature phase diagram of the 
system using real-space dynamical mean-field theory in combination with the quantum Monte Carlo 
method. We find that regardless of the intermediate dimensions examined, the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state survives to reach about one third of the BCS critical temperature of 
the spin-density balanced case. We show how the gapless nature of the state found is reflected 
in the local spectral function. While the FFLO state is found at a wide range of polarizations 
at low temperatures across the dimensional crossover, with increasing temperature we find out 
strongly dimensionality-dependent melting characteristics of shell structures related to harmonic 
confinement. Moreover, we show that intermediate dimension can help to stabilize an extremely 
uniform finite temperature FFLO state despite the presence of harmonic confinement. 
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Fermion pairing in the presence of spin-density im- 
balance has been a fundamental issue in many strongly 
correlated systems of many fields ranging from super- 
conductors to ultracold atomic gases and neutron stars 
[![. While a large magnetic field is detrimental to BCS 
superconductivity 0, Q, it has been predicted that a 
more exotic pairing mechanism would maintain Cooper 
pairs coexisting with finite spin-density imbalance. The 
Fulde-Ferrcll-Larkin-Ovchinnikov (FFLO) state suggests 
a promising scenario where, with spin polarization, 
Cooper pairs carry nonzero center-of-mass momentum, 
exhibiting a spatially oscillating order parameter 0, Q ■ 
In three-dimensional continuum, the previous mean-field 
studies predicted only a tiny FFLO area in the phase 
diagram @. Indeed, the FFLO state remains elusive in 
experiments in spite of indirect evidence of its existence 
reported in several fields [7- 11 1. However, the stability of 
the FFLO state has been suggested to depend rather sen- 
sitively on system settings, such as the presence of lattices 
[12-15 1 an d the dimensionality T^- 22 1 . In particular, it 
has been anticipated that the FFLO signature would be 
much more visible in a dimensional crossover regime be- 
tween one-dimensional (ID) and three-dimensional (3D) 
systems where the strong 1D-FFLO character at zero 
temperature could be further stabilized by the long-range 
order supported by higher dimensions. 

In this Letter, we provide finite temperature phase di- 
agrams of spin-polarized Fermi gases in lattices of inter- 
mediate dimensions in a 1D-3D crossover regime. In- 
termediate dimensions are accessible in ultracold atomic 
gases by controllin g op tical lattices as realized for weakly 
coupled ID tubes (ll| and chains [23|. Previous works 
on the FFLO state in the dimensional crossover regime 
are done at the mean-field level for coupled tubes [l9| 
or at zero temperature for two-leg ladders 24[ and lat- 



tices [25j . However, finite temperature effects are of fun- 
damental importance and need to be fully understood 
to form a precise picture on the observability of exotic 
paired states, especially regarding shell structures of dif- 
ferent phases that occur because of the overall harmonic 
trap present in ultracold gas experiments. While at close 
to ID, the scaling of critical temperature was studied 
by the effective field theory [joj . to our knowledge the 
finite temperature phase diagram in the intermediate di- 
mension has not been systematically approached so far. 
By employing a real-space variant of dynamical mean- 



field theory (DMFT) |22|, l2a l2J with continuous-time 
auxiliary- field quantum Monte Carlo method [28l - l30| , we 
explicitly consider the effects of local quantum fluctua- 
tions at finite temperatures beyond the mean-field level, 
and the presence of a trap potential which is essential to 
the shell structures. 

We consider a trapped, attractively interacting two- 
component Fermi gas in an optical lattice of ID chains 
which are coupled to form an anisotropic cubic lattice 
(see the inset of Fig. [Ha)). For deep lattice potentials, 
this system is described by the Hubbard Hamiltonian 

ft =-*|| ^2(4la C (z+l)la + h -C-) - t± C L C «V 
ila (W) icr 

il ila 

Here a =t,i denotes the (pseudo) spin state, while I 
and i stand for the chain and the lattice site within the 
chain, respectively. The fermionic annihilation and cre- 
ation operators are c ilcr and ct , and hu± is the density 
operator. We consider a harmonic potential V* = w 2 i 2 /2 
along the chains. The model is parameterized by the 
interaction strength U, the trap frequency u)\\, the spin- 
dependent chemical potential \i a and the intra- and in- 
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FIG. 1. The phase diagram of the spin-polarized Fermi gas for different dimensionalities with two representative density and 
order parameter profiles along the trapped axis for each value of t±. In the phase diagrams N stands for normal state and pSF 
for polarized superfluid while cN refers to a shell structure where the system is in the normal state in the middle of the trap with 
polarized superfluid on the edges. Similarly, cFFLO refers to the FFLO state in the middle of the trap and polarized superfluid 
on the edges. Each errorbar is determined by the closest well-converged simulation on each side of the phase boundary while 
the boundary itself is given by the mean of these two points, (a) The phase diagram for quasi-lD lattice with t± = 0.2 with 
(b) the FFLO state and the (c) the cFFLO state, (d) The phase diagram for an intermediate interchain hopping t± = 0.4 with 
(e) the FFLO state melting to (f) polarized superfluid phase at constant polarization, (g) The phase diagram for a quasi-3D 
geometry with t± = 0.8. Panels (h) and (i) demonstrate how the FFLO state is affected by the increasing temperature. The 
inset of panel (a) is a schematic of the system geometry. All energies and temperatures are in units of i n = 1. 



terchain hoppings t|| and respectively. All energies 
and temperatures are in units of and we set i|| = 1. 
Varying the interchain hopping t± from zero to one the 
system undergoes a dimensional crossover from a collec- 
tion of ID chains to a 3D (cubic) lattice. In the calcula- 
tions, the trap frequency is set to uin = 1.1 x 10~ 2 . The 
chemical potentials are varied to control the polarization 
P = (Nf- iV^)/(A r t + N±) while keeping the total par- 
ticle number approximately constant. The interaction 
strength U is chosen for each value of t± to correspond 
to the lattice equivalent of the unitarity limit |3lj . and 
thus, the interaction strength is fixed by a fundamental 
two-body property. A further investigation of an optimal 
interaction strength for the realization of the FFLO state 
remains beyond the present work. 

In DMFT, the self-energy of the system is taken as site 



diagonal, i.e. fan) = $l,l''£i,lfan) ■ The sys- 

tem is considered to be homogeneous in the interchain 
direction, and therefore, the self-energy becomes inde- 
pendent of the chain index I. Thus, the Green's function 
of the system can be written as 

[G~ 1 (k±;iu n )]ij = [G^fan)]^ 1 - [e k± cr 3 + Y^ifa n )]dij 

in which Gjj is the non-interacting Green's function of 
a single chain, ui n is the Matsubara frequency and a is 
the Pauli matrix. The transverse kinetic term is given by 
the dispersion e k± = — 2tj_ (cos fc^+cos k y ) with the trans- 
verse quasi-momentum k± = (k x ,k y ). In this notation, 
the bath Green's function of the DMFT calculations is 
given as [Q^fan)]^ 1 = [J2k ± Guik^iujn)]^ 1 + 'E i fa n ). 
The pairing order is considered within the Nambu for- 
malism, and the order parameter is defined as Aj = 



3 



We present the phase diagrams for a quasi- ID sys- 
tem with ij_ = 0.2, a system of intermediate dimen- 
sionality with tj_ = 0.4, and a quasi-3D system with 
t± = 0.8 in Fig. [T] We find that throughout the di- 
mensional crossover, the ratio of the maximum FFLO 
critical temperature and the balanced BCS critical tem- 
perature is 7Fp™^/Tg Cg « 1/3. Taking the temperature 
of 0.7 Tfflo as a reference point, we find that the po- 
larization window for the FFLO phase grows gradually 
towards the quasi-3D limit from a value of SP = 0.06 at 
t± = 0.2 to SP = 0.10 at t± = 0.8. 

In the quasi-lD regime we find a supcrfluid order pa- 
rameter which has its maximum value away from the cen- 
ter of the trap; this is clearly visible in the FFLO state 
of Fig. [lib). in this regime the FFLO state melts to 
the shell structure displayed in Fig. [TJc) in which there 
is a polarized superfluid on the edges of the trap and 
an oscillating order parameter at the center, labeled as 
cFFLO in Fig [Ha). However, above P = 0.18 the sys- 
tem reaches the normal state in the middle of the trap 
(cN) before the cFFLO shell structure is visible. The 
maximum critical temperature of the cFFLO phase is 
^cFflo = 0-45 ?b CS . Below the polarization of P = 0.13 
the cFFLO shell structure melts further to a polarized 
superfluid phase, and above P = 0.13 to the cN shell 
structure with the normal state in the middle of the trap. 

The phase diagrams for systems of intermediate di- 
mensionality (t± = 0.4) and for quasi-3D (t± = 0.8) are 
qualitatively similar. In particular, we always find the 
strongest pairing in the middle of the trap similar to a 
3D system. Comparing Figs. [H[e) and QJf ) we see how, 
here in the case of t± = 0.4, the FFLO state melts to 
the polarized superfluid phase at constant polarization. 
It is noteworthy, that in the pSF phase the edge of the 
supcrfluid still exhibits an FFLO-type oscillating order 
parameter, which gradually disappears with increasing 
temperature. As indicated in Figs. (Ud) andQJg), there 
is a transition from FFLO to the normal phase above 
the polarizations P = 0.27 and P = 0.34 for t± = 0.4 
and t± = 0.8, respectively. Figs. [T^h) and HJi) demon- 
strate how the FFLO state responds to an increase in 
temperature at constant polarization; the paired region 
recedes towards the middle of the trap. Furthermore, 
the magnitude of the order parameter decreases while the 
wavelength of the FFLO oscillation grows; this observa- 
tion holds at all dimensionalities. At the low temperature 
limit of the phase diagrams, we find good agreement with 
the zero temperature results of [25| . 

Throughout the crossover, the lower boundary of the 
FFLO region in the phase diagram with respect to po- 
larization increases with temperature. This behavior is 
essentially explained by the fact that at higher temper- 
atures, the polarized superfluid accommodates a larger 
spin-density imbalance within its thermal excitations. In 
the shell structure of a trapped gas this leads to an ad- 
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FIG. 2. At the crossing point between quasi-lD and 3D-like 
behavior, at tj_ = 0.3, we find that the FFLO state (here with 
P — 0.17 and T — 0.04) exhibits an exceptionally uniform 
oscillation amplitude across the whole system in spite of a 
density profile varying with the trap potential. 



ditional effect which counteracts to the FFLO instability. 
Namely, the polarization can be redistributed towards the 
BCS-like regions with increasing temperature. On the 
other hand, higher temperatures are less favorable for any 
pairing effects to take place and thus the upper bound- 
ary of the FFLO region with respect to polarization is 
diminished by temperature. Consequently, throughout 
the crossover the highest attainable critical temperature 
for FFLO is rather sensitive to polarization. 

Intriguingly, we find that the amplitude of the FFLO 
order parameter is essentially uniform, as shown in Fig. [21 
in a broad polarization and temperature range at the in- 
terchain hopping of t± = 0.3. This is also the crossing 
point between ID-like and 3D-like physics. For instance 
at the temperature of T = 0.05, such uniform behav- 
ior occurs in a 60 % interval of the FFLO polarization 
range in the phase diagram which is at this temperature 
5P = 0.09. This suggests that the "sweet spot" antici- 
pated in [l9| for observing the textbook FFLO order pa- 
rameter resides at tj_ « 0.3. The uniformity can also be 
an advantage for observing the state, considering probes 
that rely on strict periodicity of the order parameter. 

Further characteristics of the FFLO state can be in- 
ferred from the local spectral function plotted in Fig. [3] 
The local spectral function is defined as A, jCT (u;) 



-21mG jj<7 (iu! n 



i0 + ), and can be interpreted as 
the local density of states. We use the maximum entropy 
method to carry out the analytical continuation from the 
on-site Green's function obtained from the QMC solver 
32| . From the spectral function one can clearly see that 
the FFLO state is gapless, and also in this sense the well- 
known mean-field characterization of the state remains 
valid. Moreover, we find that only the energy states very 
close to the Fermi level contribute to the formation of 
the FFLO density oscillation. 

In conclusion, our work, which incorporates both fi- 
nite temperature effects and local quantum fluctuations, 
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FIG. 3. (a) The FFLO state near the trap center with 
t ± = 0.8, T = 0.067 and P = 0.30. The local spectral func- 
tions of the FFLO state on the lattice sites b and c of panel 
(a) for a 4- are plotted in panels (b) and (c), respectively. 
The oscillation of the density profile is revealed to affect the 
most the particles at energies close to the Fermi surface. The 
gapless nature of the FFLO state is clearly visible in the ma- 
jority (t) component on site c i.e. at the node of the order 
parameter. 



shows that the FFLO state is significant throughout a 
dimensional crossover between ID and 3D lattices at fi- 
nite temperatures. The critical temperature of the FFLO 
state is approximately one third of the superfluid critical 
temperature regardless of the dimensionality, reaching 
values as high as T ~ 0.13 in. We find that dimensional- 
ity has a clear effect on the melting behavior of the shell 
structures in a trap, which is essential in distinguishing 
between the different phases. Furthermore, we identify 
t± = 0.3 as the dimensionality crossover point that pro- 
vides the sweet spot of observing a uniform FFLO order 
parameter despite the harmonic trap confinement. On a 
final note, our results confirm that, even at the presence 
of local quantum fluctuations, the FFLO state has a wide 
region of stability in a lattice, which is in stark contrast 
to the theoretical predictions in free space where the pa- 
rameter area for the FFLO state is vanishingly small. 
This gives evidence to the fundamental role the Fermi 
surface shape has in stabilizing exotic superfluidity, and 
brings about a significant degree of freedom to future ex- 
periments aimed to realize elusive phases of matters such 
as the FFLO state. 
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